Measurement of energy eigenstates by a slow detector 
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We propose a method for a weak continuous measurement of the energy eigenstates of a fast 
quantum system by means of a "slow" detector. Such a detector is only sensitive to slowly-changing 
variables, e. g. energy, while its back-action can be limited solely to decoherence of the eigenstate 
superpositions. We apply this scheme to the problem of detection of quantum jumps between energy 
eigenstates in a harmonic oscillator. 



Observation of quantum mechanical behavior of 
macro-scale systems is one of the most fascinating cur- 
rent problems in physics. Several successes have been 
claimed so far. They include observation of avoided 
crossing between counter-propagating current states in 
Josephson loops coherent evolution of macroscopic 
quantum states in charge 0] and phase Q qubits, and 
interference experiments on very large molecules 

Recently, experiments on fabricated nano-mechanical 
systems became feasible |5[, with the position measure- 
ment resolution approaching the zero-point motion un- 
certainty Q. However, the commonly used linear posi- 
tion measurement has fundamental limitations due to the 
measurement apparatus back-action, which ultimately 
drives the oscillator away from the quantum regime 
@; B Q • Thus a question naturally arises whether it is 
feasible to observe quantum behavior of nanomechanical 
systems by means of a non-linear measurement proce- 
dure. Specifically, since the oscillator energy is quadratic 
in displacement x, coupling the measurement apparatus 
to x 2 can naturally lead to the measurement and projec- 
tion into the energy eigenstates. The main difficulty with 
this approach is that such coupling requires very precise 
tuning f]~oj ] in order to eliminate the linear-x coupling 
which would otherwise dominate at small displacements. 

There is an intimate connection between the measure- 
ment process and interaction with a general extrinsic en- 
vironment since both represent coupling of a system to 
a macroscopically large number of degrees of freedom. 
The practical difference is that the classical output of the 
measurement apparatus (or "meter") is available to the 
observer, and thus one can ask questions, e.g. about the 
signal to noise ratio, that is how much of the output can 
be attributed to the system that one is measuring, and 
how much to the noise of the apparatus Depend- 
ing on the relative strengths of the system's coupling to 
the environment and to the measurement apparatus, its 
own behavior can be either predominantly affected by 
the environment (eins election) [l2| or by the meter. The 
latter case is typically equivalent to coupling to a non- 
equilibrium environment, which nevertheless can result 
in a crossover to effectively thermal behavior of the sys- 
tem with the "temperature" determined by the degree 
of non-equilibrium in the meter (e.g. bias voltage in a 



quantum point contact [8|]). 

When environment is weakly coupled to the system 
and is slow compared to the intrinsic system dynamics 
(determined by its self-hamiltonian Hq), an interesting 
possibility arises: While the environment can lead to the 
decoherence of superpositions of the energy eigenstates, it 
can not effectively cause transitions between them. This 
is a consequence of the quantum adiabatic theorem. De- 
coherence is commonly associated with a gradual projec- 
tion into one of the states of the superposition. Thus one 
concludes that the system should remain in the energy 
eigenstates [13| most of the time, with the rare transi- 
tions (quantum jumps, QJ) between them governed by 
the residual relaxation processes. A necessary condition 
for this scenario is a presence of non-zero diagonal (in the 
energy eigen-basis) coupling between environment and 
the system, which leads to environment-induced fluctua- 
tions in the system level spacing. 

In this Letter we apply the ideas of slow environment 
to the measurement process and determine the circum- 
stances under which QJ should become observable. We 
perform analysis for a harmonic oscillator playing the 
role of the quantum system being measured, with the 
Hamiltonian 



P 
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where x is the oscillator position, p is the momentum, 
to is its mass and luq is its frequency. Analysis for a 
more general systems is similar. For harmonic oscilla- 
tor, the requirement on non-zero diagonal coupling is 
equivalent to inclusion of non-linear terms in the meter- 
oscillator coupling. For concreteness we consider terms 
up to quadratic order in displacement, 



-Hint = (Aix + \ 2 x )F = Sl{x)T. 
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Note that for a slow operator T the linear coupling term 
(Ai) does not generate non-linear in x terms within per- 
turbation theory. 

We will first derive the measurability criteria for QJ 
for a general measurement apparatus that interacts with 
the oscillator via operator T and generates output /. We 
will then give a specific example of such a measurement 
when a single electron transistor (SET) is capacitively 
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coupled to the oscillator. This experimental procedure 
is closely related to the displacement measurement tech- 
niques currently in use Q . 

The complete Hamiltonian for the combined system 
of the oscillator and the measurement apparatus (meter) 
consists of three parts, 



H = H„ 
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with H a and H lnt specified above and Hm being the 
Hamiltonian of the meter. The evolution of the full den- 
sity matrix p is given by p — i[p, H]. Assuming that the 
interaction between the measurement apparatus and the 
system is relatively weak, that is it dominates neither 
the oscillator's dynamics nor the meter's ("weak mea- 
surement"), we can determine the evolution by means 
of the perturbation theory. In the interaction represen- 
tation, where p = e -K^o+H M )t pe i(H +H M )t anc j jj. nt = 
e i(.H +H M )t H .^ e -i(H +H M )t^ the secon d-order perturba- 
tion theory yields, 

'p = i\p,H iB/t ]=i\p(0) + i f dt%H- mt } t ,,H- mt ] (3) 

Jo 

« i[p(0),H int \- f dt'[[p,H int (t%H int ] (4) 
Jo 

Here the operators with unspecified time variable are 
taken at time t. The final step is equivalent to the Marko- 
vian approximation, which is valid if the time correla- 
tions in Hi n t decay faster than the decay times of di- 
agonal and off-diagonal elements of the density matrix 
(in the jargon of nuclear magnetic resonance, the times 
T\ and T2, respectively). Making now the common as- 
sumption that the initial density matrix is factorizablc, 
p(0) = p o (0) 0Pm(O), and tracing over the meter den- 
sity matrix pm we can obtain the evolution of the reduced 
density matrix of the oscillator p a . 

' Po = i[p ,n](t) (s) 

rt dt'[[p .n(t%n]^^l 
I. & 

1 dt'[{p ,Ci(t')}A^^ 

The first term in Eq. ^ corresponds to the renormal- 
ization of the oscillator Hamiltonian by the interaction 
with the meter via the average (.F(i)) = TtmPmF^S) and 
will be neglected in the following discussion based on the 
assumption that the interaction between the meter and 
the oscillator is week. The second and the third terms in 
Eq. ([5]) describe decoherence and relaxation of the oscilla- 
tor induced by the backaction of the meter. Defining the 
noise correlator G F (t,t') = {T{t)T{t')) - (T{t))(F(t')), 
the symmetric and antisymmetric correlators S F and A F 
are G F (t, t')±Gp(t', t). Note that the antisymmetric cor- 
relator is only non-zero if the operator T is quantum, i.e. 
doesn't commute with itself at different times. 



The density matrix evolution equation can be simpli- 
fied by substituting fi(i') = e lH " {t '- t) n(t)e- tH ^ t '" t '> . As- 
suming the stationary noise correlator G(t, t') = G(t — t') 
which decays on the time scale faster than T\ and T 2 , in 
the energy eigen-basis, 
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hnn = -A^[^(o; )(n + l/2)-i^(a;o) 
+\\tlGjr(uj )(n + 1) p n+1>n+1 
+X 2 £ 2 l G?-(-u!o)n p n -i, n -i 

where n ^ fc in the first equation, and tilde denotes the 
Fourier transform. The zero-point motion amplitude is 
£0 = ■wh/ {2midQ). From Eq. ([7]) the relaxation rate from 
state n is 
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We can similarly define from Eq. ([6]) the decoherence rate 
corresponding to the decay of the off-diagonal density 
matrix element p nm , 
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For slow measurement we expect that Tmk <C Ti„ for 
any n 7^ fc. 

For the energy eigenstate n to be measurable, the re- 
laxation time for that state, T\ n , has to be longer than 
the measurement time Tm- Typically Tm is defined as 
the time needed to measure a variable n with the signal- 
to-noise ratio of 1. For low- frequency part of the detector 
output noise 5/(0), and the incremental change in the 
output AI = \I(n + 1) — I(n)\, the measurement time is 



Tm = 



5/(0) 

(Air 
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From the linear response theory, we can define the 
forward A = (i/K) J™ dt[T(Q), I(t)] and reverse A' = 
(i/h) f™dt[I(0), F(t)] detector gains. Then 



AI = \ 2 ilkAn = A 2 ^A. 



(11) 



The signal-to-noise ratio (SNR) achievable before the n- 
th state decays is 
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To derive the last inequality we used the Schwartz 
inequality, SiSy < \Sir\ 2 , the identity A - A' = 
(2/H) Im5/jf, and the absence of positive feedback in the 
detector, A • A' < 9]. The conditions for the equality 
in Eq. (|13|) include absence of reverse gain, A' = 0, and 
ReSir = 0. For realistic nanomechanical systems' fre- 
quencies and temperatures, the first factor can be of order 
1 , the second is likely to be small due to the square of the 
zero-point-motion amplitude, and the last one, by con- 
struction, for the slow detector can be very large. Thus 
as a matter of principle the measurement of the number 
states by a slow detector is possible. 

In addition to the coupling to the meter, an oscillator 
is influenced by its own environment, which leads to a 
finite intrinsic quality factor Q. The contribution of this 
heat bath to relaxation at T > Hujq is 



1 u 



T 



(14) 



which may reduce the SNR of Eq. (JT3 
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FIG. 1: Simulation of the oscillator occupation number dy- 
namics, described by Eq. (0. The parameters of simulation 



are r up — 1, r^ n — 1.2 



The output of a detector with high SNR, r\ n 1, is 
proportional to nit) (QJ's are not masked by the de- 
tector noise). Stochastic simulation of this process is 
shown in the top panel of Fig. [1] The switching be- 
tween the plateaus corresponding to different values of 
n occurs with the average rate F ~ (T up + T c i n )n — 
r up (r up + T dn )/{Tdn - r up ), where T up = A^oG^Wo) 
is the transition rate from the ground state into the 
first excited (0 — * 1), and Td n = A^qGjf(— uiq) is the 
opposite rate (1 — * 0), sec Eq. (|7|). The fast switch- 
ing is superimposed on top of the slow dynamics occur- 
ring on the time scale of the inverse damping coefficient, 

7 _1 = (Tdn - Tup)" 1 . 

In the presence of a detector noise the QJ may not 
be directly observable in I(t). However, there still 



may be signatures in the signal correlations. First 
let's consider the two-point time autocorrelation func- 
tion C 2 (ti,f 2 ) = (n(ti)n(t 2 )) which enters the output 
autocorrelation function, (I(ti)I(t2)). It can be ex- 
pressed in terms of the stationary probability distribu- 
tion P(n) and the conditional probability that the sys- 
tem is in state n 2 at time t 2 given that it was in state 
n\ at an earlier time t\, P(n 2 , i 2 ; ni, ii), as C 2 (ti,i 2 ) = 
Zmi,n2 n 2P( n 2,h;ni,ti)niP(ni). From the equation of 
motion for n(t) = Y] np nn (t) that can be obtained from 
Eq. 0, it is easy to see, however, that C2 is not sensitive 
to the fast switching with rate f . 




FIG. 2: Four-point correlation function containing the infor- 
mation about fast switching between the number states. The 
parameters of simulation are T up = 1, Tdn = 1.2 

Therefore, in order to extract the signatures of the fast 
switching, we need to go to higher-order correlators. One 
approach is to convert the multi-level staircase dynamics 
into a two-level, or telegraph switching. This is achieved 
by performing the following non-linear transformation on 
the output signal: I(t) -> [I(t + r) - I{t)] 2 . It will be 
telegraph-like if Tt <C 1. Examples for different values of 
r are shown in the lower three panels of Fig. [TJ We now 
define the four-point correlation function, 

C 4 (ti,t 2 ;r) = ([n(t2+r)-n(t 2 )] 2 Kti+r)-n(ti)] 2 ), 

which is expected to be sensitive to the fast switching. 
The contributions to this correlation have the general 
form {n{ti)n{t^)n{t2)n{ti)) , with t± > t 3 > t 2 > t\- Sim- 
ilar to the two-point correlation function, it can be ex- 
pressed as 



E 



n 4 P(ri4, t 4 ; n 3 , t 3 )n 3 P(n 3 , t 3 ; n 2 ,t 2 ) (15) 
xn 2 P(n 2 , i 2 ; ni,ii)niP(ni). 



Rewriting Eq. ([7]) in the matrix form, p — Mp, the tran- 
sition matrices can be expressed as 

P(n2,t 2 ;ni,ti) = Uexp[D(h - ti)p-\ (16) 
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where D is a diagonal matrix of eigenvalues and U is the 
matrix of eigenvectors of M, MU — UD. Performing nu- 
merical diagonalization on the truncated system of equa- 
tions ([7]), we obtain the correlator C 4 (t;r) = C±(t,0; r) 
for various values of r, as shown in Fig.O Indeed we find 
a qualitative change when the delay time r becomes com- 
parable to the average switching time, Tt f» 1. For short 
delay times, Tt 1, the correlation function rapidly de- 
cays at short times (t < r), followed by a slow exponential 
decay at long times. It is easy to see that in this regime, 
the correlation function has to be linear, C4 oc r — t for 
t < t. For long delay times, Tt > 1, the correlation func- 
tion is a sum of exponentials. When an additive detector 
noise £ with short correlation time << F -1 is included, 
the behavior of C4 only changes at t < t£. Thus, as far 
as T£ < r, the extrinsic noise can be effectively filtered 
out and C4 can be used to detect QJ's. 
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FIG. 3: Scheme for SET displacement detector. Position of 
the oscillator modulates capacitance C g (x), thus modifying 
the current between the left and right leads. 



SET is a nearly ideal detector |l4|; thus, the conditions 
for equality in Eq. (|13p are approximately satisfied. As- 
suming that the intrinsic (measurement-unrelated) relax- 
ation of the oscillator is smaller than the measurement- 
induced one, the signal-to-noise ratio is 
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From Eqs. ([S]) and (|14p the measurement-induced 
relaxation will dominate the environmental one if 
\{llQ/{TT) > 1. 

For a numerical estimate we take an oscillator with 
the frequency ojq = 2tt x 100 MHz and the mass m — 
2.5 • 10~ 17 (for 1 /im x 0.1 /.tm x 0.1 /im dimensions at 
the Si density 2,500 kg/m 3 ). The zero-point-motion am- 
plitude is £q = 10~ 13 m. Assuming that the effective 
separation between the nanoresonator and the SET is- 
land is d = 1 nm and T = 10 4 s _1 , for SNR we find 
Vn ET ~ 100/(2n+l). It is easy to check that the intrinsic 
relaxation (finite Q) under these conditions is insignifi- 
cant compared to the SET-induced relaxation. At tem- 
perature T = 10 mK (superconducting SET) the thermal 
occupation number of the resonator is about 2. There- 
fore, under these conditions one can expect to observe 
QJ between values of current that correspond to differ- 
ent energy eigenstates of the nanoresonator. 
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A particular example of a measurement apparatus that 
can be used to perform the slow measurement described 
above is a Single Electron Transistor (SET). The nanores- 
onator interacts with the SET island via displacement- 
dependent gate capacitance C g {x) = Cg/(l — x/d), where 
d is the distance between the nanoresonator and the SET 
island (Fig. [3]). The displacement is read out by measur- 
ing the current / through the SET. The island is also 
coupled to the leads via capacitances C\ and C2. For 
concreteness assuming that C° = C\ + C2, the interac- 
tion Hamiltonian is 

Ws + aO' < 17 > 

which indeed has the form of Eq. ([1}, with T = q the 
island charge, Ai = V g /(4d) and A 2 = V g /{8d 2 ). To 
satisfy the conditions of "slow" measurement, the SET 
tunneling rate has to be much slower than the oscillator 
frequency, r <C u>o- We consider the "threshold" regime, 
which corresponds to the resonant level nearly aligned 
with one of the leads' chemical potential and corresponds 
to the highest read-out sensitivity. We assume that cither 
T <C T or that SET is superconducting. Then, 

S q (0)~L S q (u )~^. (18) 
1 iv 
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